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We study the long time coherence dynamics of a two-mode Bose-Hubbard model in the Josephson 
interaction regime, as a function of the relative phase and occupation imbalance of an arbitrary 
coherent preparation. We find that the variance of the long time fluctuations of the one-body 
coherence can be factorized as a product of the inverse participation number 1/M that depends only 
on the preparation, and a semi-classical function C(E) that reflects the phase space characteristics 
of the pertinent observable. Temporal fluctuations can thus be used as a sensitive probe for phase 
space tomography of quantum many-body states. 
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The two mode Bose-Hubbard Hamiltonian (BHH) ap- 
pears in different guises in a perplexing variety of fields. 
Cast in spin form, it is known in nuclear physics as the 
Lipkin-Meshkov-Glick (LMG) model of shape phase tran- 
sitions [H- It is broadly used to describe interacting 
spin systems and serves as a paradigm for squeezing 
and entanglement Q. As such, it offers schemes for the 
generation of squeezed states for optical interferometry 
below the standard quantum limit and its matter- 
wave equivalent . It is commonly employed to describe 
the Josephson dynamics in systems of bosonic atoms in 
double-well potentials @ and suggests prospects for the 
generation of macroscopic superposition states @- The 
same model is also known in condensed matter physics 
as the integrable dimer model @ with applications to 
the dynamics of small molecules, molecular crystals, and 
self-trapping in amorphous semiconductors. 

Like the paradigmatic Jaynes-Cummings model in 
quantum optics joj, the bimodal BHH dynamics with a 
coherent spin state preparation exhibits a series of col- 
lapses and revivals of its single-particle coherence due to 
interactions [lol - [l3j . These recurrences are manifested in 
the collapse and revival of the Rabi- Josephson popula- 
tion oscillations, or of the multi-realization fringe visibil- 
ity, when the two condensates are released and allowed 
to interfere. Below we study the long time BHH dy- 
namics for general coherent spin preparations \0,<p). In 
such states all particles occupy a single superposition of 
the two modes, with a normalized population imbalance 
S z = cos(#) and a relative phase <f>. 

The characteristics of this dephasing-rephasing dy- 
namics strongly depend on the dimensionlcss interac- 
tion parameter u = UN/K, where U is the interaction 
strength, TV is the total particle number, and K is the 
hopping amplitude. In the linear Rabi regime (u < 1) 
time evolution is straightforward because the interaction 
is not sufficiently strong to modify the one-body nature 
of the dynamics. Accordingly, one observes only coher- 
ent Rabi oscillations in the population difference with a 
typical frequency cjj = ^/K(K — UN) = Ky/1 — u m K 
which reflects mainly the coupling K between the two 
modes, accompanied by a small modulation of the single- 
particle coherence at the same frequency. 

The coherence dynamics in the strong interaction Fock 
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FIG. 1: (color online) Collapse and revival dynamics of 
the normalized population imbalance S z (top) and of the 
single-particle coherence S (bottom) for representative equal- 
population coherent preparations \tt/2,(/>), with cj> — (o), 
0.2% (*), 0.4tt (V), 0.6vr (o), 0.8?r (□), and % (A). Note that 
in the last case S z (t) — at any moment, while S(t) oscillates 
wildly. The BHH parameters here and in all subsequent fig- 
ures are N = 100 and u = 5, corresponding to the Josephson 
regime. 



regime (u > N 2 ) is also fairly simple because it reflects 
the Fock basis expansion of the initial coherent prepa- 
ration. The expected coherence recurrences have been 
observed experimentally [3, EH with a striking demon- 
stration of exceptionally long time dynamics, allowing to 
probe effective multi-body interactions through the de- 
pendence of U on the number of atoms [l6[ . 

The dynamics in the Josephson regime (1 < u < N 2 ) is 
by far richer and more intricate, reflec ting the coexistence 
of three distinct phase space regions [ill. Il2j. Unlike the 
Fock-space recurrences which only depend on the pop- 
ulation imbalance, the Josephson coherence dynamics is 
also highly sensitive to the relative phase. Previous work 
has considered specific preparations that were of contem- 
porary experimental relevance, e.g. small perturbation of 
the ground state that results in Josephson oscillations, or 
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FIG. 2: (color online) Long time average (a,b) and variance 
(c,d) of the one-particle coherence S(t), starting from the co- 
herent preparation \6,cj>). Panels (b) and (d) show a section 
along the Equator 9 — tt/2 with <\> € [0, ir] (blue solid line), 
with appended sections along the arcs <j> = with 6\ G [tt/2, 0] 
(black dashed line), and <f> = n with 62 £ [0, 7r/2] (red dashed 
line). Symbols correspond to values which were calculated 
from a direct numerical propagation of the states in Fig[Tj 

a large population imbalance that leads to self-trapping 
[17( 1 . Here we adopt a global, tomographic approach by 
characterizing the long time temporal quantum fluctu- 
ations for all possible coherent preparations. This ap- 
pears to be a formidable task, but as shown below, a 
relatively simple semi-classical perspective provides an 
adequate framework for the required analysis. 

Observations.— The intricacy of the Josephson 
regime quantum dynamics is illustrated in FigJT] and in 
FigEJ Various coherent preparations lead to qualitatively 
different behavior, depending on the initial population 
imbalance cos 9 and on the relative phase <fi. The cause 
of this diversity is that different coherent preparations 
sample different parts of the spectrum. Each coherent 
spin state constitutes a superposition of eigenstates that 
can be associated with qualitatively different regions in 
the corresponding classical phase space. This is in a stark 
contrast to the Rabi and Fock regimes, where the eigen- 
states occupy, so to say, a single component phase space 
that allows only one type of motion. 

Our main finding below is related to a striking fac- 
torization in the calculation of the variance of the tem- 
poral fluctuations of the pertinent observables. Given 
an observable A with a fluctuating expectation value 
A{t) = (A) t that has long time average A(t), we find 
that its variance can be approximated quite accurately as 
a product of a quantum term and a semi-classical term: 

a] e If-I(i) 2 = Ji CA{E) - (1) 

Here M is the participation number which depends only 
on the preparation, while Ca (E) corresponds to the clas- 



sical fluctuations of A along a mean-field trajectory that 
has an energy E. The implications of this observation 
for phase space tomography arc further discussed in the 
concluding paragraph. We now turn to a detailed pre- 
sentation of the analysis. 

The BHH. We consider a similar scenario to that 
in Rcf. [15(, which observed the long time collapse and 
revival of coherence in the Fock regime. For such strong 
interactions the modes are essentially decoupled, hence 
the dynamics is fully captured by the Gutzwiller ansatz 
of a direct product of single-site states, each of which is 
a coherent wavepacket of number states [H], EE EE EH ■ 
Thus, in the Fock regime the two-mode BHH generates 
precisely the same coherence dynamics as the many-site 
BHH of a BEC in an optical lattice. This allows for mon- 
itoring the fringe visibility in single shot intcrfcromctery 
of an optical lattice, rather than repeating a two-mode 
experiment many times [HI, EE EH- In the Josephson 
regime, the dynamics of a lattice with one mode per site 
is quite different from the two-mode dynamics. However, 
the two-mode model can be realized with two spin com- 
ponents in each isolated site [2fj| or with an array of inde- 
pendent double wells [2l| , thus retaining the convenience 
of single-shot measurements. 

Assuming that no bias field is applied, the pertinent 
BHH is, 

% = — K (h\a 2 + a^a-ij 

+y [hi {hi - 1) + h 2 (n 2 - 1)] , (2) 

where hi and a\ are bosonic annihilation and creation 
operators, respectively. The particle number operator in 
mode i is fii = a\a,i. Since the total particle number 
h\ + h% = N is conserved, we can eliminate respective 
c-number terms and obtain the BHH in spin form, 

H = -KJ X + UJ 2 Z , (3) 

where J x ={a\a2+a\a\) /2, J y ~(a\a2~ a\ai)/(2i), and 
Jz=(hi—n2)/2. The number conservation becomes an- 
gular momentum conservation with j = N/2. Below we 
assume for simplicity that the interaction is repulsive 
U > 0, but the U < case amounts to a simple transfor- 
mation K H> —K, and E i-> —E. Thus the phase space 
with attractive interaction is simply an inverted mirror 
image of the repulsive-inetraction case and there is no loss 
of generality. In this spin representation, each state is 
characterized by the normalized Bloch vector S = (J) /j, 
where its z projection S z = cos(#) corresponds to the 
normalized population imbalance, its azimuthal angle 4> 
corresponds to the relative phase between the modes, and 
its length S corresponds to the single-particle coherence. 

The phase space structure of the BHH is set by the 
previously defined dimensionless interaction parameter 
u. Its characteristics in the three interaction regimes are 
discussed in great detail elsewhere [lEUE!- Briefly, 111 
the Josephson regime (1 < u < N 2 ) the phase space con- 
tains two 'islands' that are separated from a 'sea' region 
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by a separatrix. The sea trajectories correspond to Rabi- 
Josephson population oscillations, whereas the island 
traj ectories correspond to self-trapped phase-oscillations 
p7| . In the Fock regime (u > N 2 ) the sea becomes too 
small to support quantum states, while in the opposite 
limit - in the Rabi regime (u < 1) - the islands disappear, 
so that only Rabi-type oscillations are feasible. 

Evolution.— We study the dynamics induced by the 
Hamiltonian of Eq.([2]), or equivalcntly by Eq.Q, starting 
from an arbitrary spin coherent state preparation, 
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\0,4>) = ^ |cos(0/2)a{ +sin(6»/2)e^a + 

= exp(-i(j)J z ) cxp(-i6»Jy) exp(i<j)J z )\j, -j},(4) 

where |vac) and \j, —j) are the vacuum states of the 
Heisenberg-Weil and SU(2) algebras, respectively. We 
focus our attention on experimentally relevant observ- 
ables such as the population imbalance S z and the single- 
particle coherence S. Note that S is the best fringe visi- 
bility one may expect to measure by proper manipulation 
of the Bloch vector, i.e. if we are allowed to perform any 
SU(2) rotation. The expected fringe visibility if inter- 
fcrometry is carried out without further manipulation is 

9i2 = (l/i)[l^| 2 + l^yl 2 ] 1 ^ 2 - F° r presentation purpose 
we focus on S, but the results for are very similar. 

The evolution of the population imbalance and of the 
one-particle coherence for several representative coher- 
ent preparations is shown in FigfT] For each preparation 
\6, 4>) we obtain the long time- average and the variance of 
S z (t) and S(t). These numerical-propagation values for 
the preparations of FigQ] are represented by the symbols 
in Fig[2f b) (d). In order to see the overall picture, we plot 
in Fig[Ua)(c) an image of S(0, 4>) and as (9, <f>) for all the 
possible coherent preparations. We note parenthetically 
that in practice, the exact values are not obtained from 
finite time simulations, but deduced from the energy ex- 
pansion as detailed in the analysis section below. As is 
exemplified in Fig.[2](b)(d), the agreement between the 
energy expansion (lines) and the numerical propagation 
(symbols) clearly validates the pertinent expressions. 

Analysis.— In order to deduce the exact time average 
of any A(t), we expand it in the energy basis as 



A(t) 



I'.fl 



c* v c^A vli exp[(E v - E^t/h], (5) 



where \E V ) are the BHH eigenstates, c v = (E v \ij]) are the 
expansion coefficients of the initial state \4>) = \9, </>), and 
A vfl = (E^\A\E^). The long-time average eliminates the 
oscillating terms, hence A(t) = ~}2 v PvA yl y, with proba- 
bilities p v = |c„| 2 , while the variance is 
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The matrix elements in Eq. ((HJ) can be evaluated semi- 
classically using the following prescription [12, [22[: a 
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FIG. 3: (color online) Calculation of C(E) for the observable 
A — J z . Panel (a) is an image of the matrix elements |A nm |, 
with color-scale in loglO units. Panel (b) displays equal energy 
contours E = (E n + E m )/2. In panel (c) the function Ca{E) 
is evaluated according to the middle (symbols) and the r.h.s. 
(lines) of Eq.(f7]). The dashed vertical line marks the separatrix 
energy E — NK/2. Panel (d) is a phase space image of 
Ca(E). The left half of the image corresponds to the lines of 
panel (c), whereas the right side corresponds to the symbols. 



classical trajectory of energy E is generated using the 
BHH mean field equations of motion, and A c i(t) is cal- 
culated; then the classical power-spectrum C^(cj) is ob- 
tained via a Fourier transform of \A c i(t) — Ad]; and fi- 
nally the result is divided by the mean level spacing g 
at that energy, providing the approximation 1^4^! 2 ~ 
C^{E v —E tl )/(2-KQ). This is a very general procedure 
which is usually applied to chaotic systems, but it ap- 
plies equally well to the integrable non-linear motion of 
the two-mode BHH. Since the motion is periodic but 
highly non-linear, the integral over C^(u>) turns out to 
be the variance of a periodic function that contains many 
harmonics. One may say that the power spectrum of 
A(t), though quantized, is semiclassically narrow but 
quantum-mechanically broad. The number of participat- 
ing eigenstates that contribute to Eq.©, is convention- 
ally evaluated as M e E^P 2 -] -1 - Assuming M>1, ap- 
proximating p v ps 1/M, and neglecting non-participating 
eigenstates, one obtains Eq.(Q} where, 



c A (E)= £ kWI 2 = I Ci{u) 



M>o 



2tt' 



(7) 



Above r = [y — (i) is the diagonal coordinate of the ma- 
trix, and it is implicit that the summation is carried out 
over a section v + \x = const such that (E v + E^)/2 ~ E. 

Numerical verification.— In Fig. [3] we illustrate 
a representative evaluation of Ca(E) for the popula- 
tion imbalance A = J z . The matrix elements (Jz) v „ are 
shown in panel Fig(3](a), confirming the assumptions of a 
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FIG. 4: (color online) The participation number of the co- 
herent state preparations, and the long time fluctuations of 
the observable A = J z . Panel (a) displays the participation 
number M(9,(f>) of the states \9, <f>). Panel (b) provides a 
cross-section of M(9, <f>) along the same contour that has been 
defined in Fig.[2] Panel (c) overlays M(8,4>)o\(9,(j)) with 
the constant energy contours in phase space. Panel (d) dis- 
plays the variance a\ as calculated from the approximation 
Eq. {J} (right half of the panel), compared with the actual 
value Eq.(|6} (left half of the panel). 



semiclassically narrow but quantum-mechanically broad 
spectrum. The equal energy contours in panel Figj3]Jb) 
demonstrate that due to the semiclassical narrowness we 
can approximate the energy shell average by a summa- 
tion over the diagonals v + /x = const. The results of the 
summation over the matrix elements and the integration 
over the classical fluctuations according to Eq. ([7]) are 
compared in panel Fig[3][c), and the final contour map of 
the power spectrum in the phase space is presented for 
both methods in panel Fig[3jd) . This comparison demon- 
strates the validity of the semiclassical approximation of 
the matrix elements. 

In order to evaluate the variance of the fluctuations, we 
need to calculate the participation number M for a gen- 
eral coherent preparation \9, tt). The result is shown in 
Fig[Ua)(b). Due to the factorization Eq.([T]), this function 
needs be calculated only once for all desired observables. 
While wc do not have a closed analytic expression for 
M{6, </>), its characteristic value and its dependence on u 
and N in the different phase space regions can be evalu- 
ated from general considerations as detailed in Ref. [12]. 
Generally speaking, the highest participation numbers 
are obtained at the top of the separatrix and scale as 
M w \> r N \og(N / u) , i.e. like the square root of N. By 
contrast, the equatorial states |7r/2,0) and \tt/2,7t) have 
participation numbers of order unity: A/(7r/2,0) ss ^fu 
and M w y/u\og(N/u), respectively. 

Fig2Jc) confirms our result that the long time temporal 
fluctuations can be factorized as the ratio of the power- 
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FIG. 5: (color online) Power spectrum Cs(E) (a) and vari- 
ance of fluctuations erf (b) for the single-particle coherence 
S. Notations are the same as in FigJ3{c) and Figf3{d). Note 
that the left half of the image is the same as Fig[2jc). 



spectrum and the participation number. Indeed, when 
we plot Ma A as a function of 8 and <j> we obtain that the 
Ma\ = const contours coincide with the constant-energy 
lines. Thus, Ma\ is only a function of E. The only ex- 
ception is in the immediate vicinity of the coherent state 
|tt/2, 7r) for which the semi-classical assumption M>1 
does not hold. 

Having ascertained the semiclassical factorization, it 
remains to verify the quantitative agreement between the 
values of the contour lines of Ma\ and the semiclassi- 
cally calculated Ca(E). In Figj4|d) we use the partici- 
pation number M and the calculated Cj z (E) to predict 
the variance of the population imbalance oscillations for 
the various preparations, and compare it to the results 
obtained by a numerical propagation, or by using Eq.©. 
The agreement is good, confirming the validity of Eq.([T]) 
and demonstrating that the long time population oscil- 
lations will only survive in the vicinity of the unstable 
equal-population <f> = 7r preparation, where the power 
spectrum is large and the participation number is small. 
Note that the classical fluctuations are large for the other 
separatrix preparations, too, but away from </> = n the 
participation number is large, and hence the quantum 
fluctuations are quenched. Similar results are shown in 
Fig|5]for the variance of the coherence fluctuations. 

Summary.— The magnitude of the long-time quan- 
tum fluctuations of an arbitrary observable A can be de- 
duced from the tractable classical dynamics and from 
the a-priori known participation number of all coherent 
preparations. Moreover, because Ca{E) only depends on 
energy, it is possible to obtain a tomographic image of 
the phase space from the observed fluctuations by plot- 
ting Ma\. The contours of this plot trace the equal 
energy lines and their values are specific to the measured 
observable. Eq. (Q]) constitutes a dramatic reduction in 
complexity and offers great insight on the dynamics of 
the two-mode BHH. It is successfully employed above to 
characterize the observed fluctuations in population im- 
balance, one-particle coherence, and fringe-visibility. In 
a larger context, these findings are related to studies of 
thermalization in isolated many-body quantum systems 
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